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Abstract 

The  problem  of  partitioning  an  assembly  of  polyhedral  objects  into 
two  subassemblies  that  can  be  separated  arises  in  assembly  planning.  We 
describe  an  algorithm  to  compute  the  set  of  all  translations  separating  two 
polyhedra  with  n  vertices  in  0{n*)  steps  and  show  that  this  is  optimal. 
Given  an  assembly  of  k  polyhedra  with  a  total  of  n  vertices,  an  extension 
of  this  algorithm  identifies  a  valid  translation  and  removable  subassembly 
in  0{k^n* )  steps  if  one  exists.  Based  on  the  second  algorithm  a  polynomial 
time  method  for  finding  a  complete  assembly  sequence  consisting  of  single 
translations  is  derived.  An  implementation  incorporates  several  changes 
to  achieve  better  average-case  performance;  experimental  results  obtained 
for  composite  objects  consisting  of  isothetic  polyhedra  are  described. 
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The  problem  of  finding  sequences  of  motions  for  the  assembly  of  a  given  object 
consisting  of  polyhedral  parts  arises  in  assembly  planning.  This  problem  can  be 
regarded  as  a  motion  planning  problem  with  multiple  moving  objects.  In  this 
general  form,  the  problem  involves  many  degrees  of  freedom.  Since  known  meth¬ 
ods  for  motion  planning  allowing  general  motions  are  exponential  in  the  number 
of  degrees  of  freedom,  it  is  useful  to  restrict  the  type  of  motion  considered.  Here 
we  will  impose  the  following  restrictions: 

•  Each  step  in  an  assembly  sequence  concerns  two  subassemblies.  Two  sub- 
assemblies  which  have  been  joined  in  a  previous  step  are  not  moved  relative 
to  each  other  in  subsequent  steps. 

•  At  each  step  in  an  assembly  sequence,  a  single  translation  moves  the  first 
subaissembly  to  its  final  position  relative  to  the  second  subassembly. 


Figure  1;  An  assembly  in  which  no  single  part  can  be  removed 


These  restrictions  embody  practical  constraints  often  imposed  on  assembly  se¬ 
quences  by  manufacturing  processes;  complicated  assembly  motions  and  oper¬ 
ations  joining  more  than  two  subassemblies  make  assembly  more  difficult  and 
raise  manufacturing  costs. 

For  assemblies  of  rigid  parts,  an  assembly  plan  can  be  obtained  by  reversing 
a  valid  disassembly  plan.  In  this  context,  we  will  address  the  following  two 
subproblems; 

1.  Given  two  polyhedra,  compute  the  set  of  all  single  translations  separating 
these  polyhedra. 

2.  Given  an  assembly  A  of  several  polyhedra,  decide  whether  there  is  a  direc¬ 
tion  d  and  a  subassembly  5  C  ^  such  that  a  translation  along  d  separates 
5  from  the  remaining  parts  /I  \  5. 

Given  a  solution  to  the  second  problem,  we  will  show  that  complete  assembly 
sequences  for  polyhedral  parts  can  easily  be  computed. 

Two  examples  for  the  second  problem  are  shown  in  figures  1  and  2.  None 
of  the  polyhedra  in  figure  1  can  be  separated  from  the  remaining  parts  by  a 
translation  involving  a  single  object,  but  there  are  subassemblies  which  can  be 
moved  simultaneously.  In  figure  2  any  subset  of  the  cubes  Pi,. can  be 
removed  by  a  single  simultaneous  translation  from  the  remaining  objects.  This 
second  example  shows  that  the  number  of  removable  subassemblies  is  exponen¬ 
tial  in  general.  Hence,  it  is  not  practical  to  compute  all  removable  subassemblies 
explicitly. 

However,  we  show  that  it  can  be  determined  in  polynomial  time  whether  a 
removable  subassembly  exists.  Specifically,  we  describe  an  optimal  algorithm 
for  solving  the  first  problem  above.  The  method  is  then  extended  to  derive  a 
polynomial  time  algorithm  to  solve  the  second  problem.  Using  this  procedure. 


Figure  2:  An  assembly  of  cubes 

we  derive  a  polynomial  time  method  for  finding  complete  assembly  sequences. 
Finally,  an  implementation  of  the  above  algorithms  and  the  results  of  various 
assembly  planning  experiments  are  described. 


1  Related  Work 

A  survey  of  earlier  methods  for  separating  sets  in  two  and  three  dimensions  is 
given  in  [14].  In  [11]  lower  bounds  on  the  number  of  simultaneous  translations 
necessary  for  separating  objects  are  derived.  Dawson  [3]  shows  that  two  or 
more  star-shaped  objects  can  always  be  separated  by  translating  the  objects  in 
different  directions  simultaneously.  In  addition,  it  is  shown  in  [3]  that  for  some 
assemblies  of  convex  polyhedra,  no  individual  parts  are  removable  by  a  single 
translation. 

Homem  de  Mello  and  Sanderson  [6]  give  a  method  to  calculate  the  polyhedral 
convex  cone  containing  the  infinitesimal  translations  allowed  by  a  set  of  planar 
contacts  in  space.  A  polynomial-time  algorithm  to  identify  subassemblies  that 
are  connected  and  can  be  translated  a  small  distance  relative  to  the  rest  of  the 
assembly  is  described  in  [16].  Both  of  these  methods  consider  only  conta'" 
and  thus  cannot  find  collision-free  extended  translations. 

Krishnan  and  Sanderson  [7]  address  problem  1  by  mapping  the  set  of  all 
unit  translations  onto  a  two-dimensional  grid,  and  marking  grid  elements  that 
correspond  to  collisions  between  two  polyhedra.  Any  unmarked  elements  then 
represent  valid  removal  translations.  However,  this  method  is  otny  accurate  to 
the  size  of  the  grid,  and  cannot  be  used  to  find  translations  i.ivolving  contacts 
between  the  two  parts. 

Pollack,  Sharir  and  Sifrony  [12]  consider  sequences  of  translations  to  separate 
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polygons.  The  algorithm  [12]  is  limited  to  planar  assemblies  of  two  parts,  but 
is  able  to  find  separating  motions  consisting  of  several  distinct  translations. 

Toussaint  [15]  describes  an  algorithm  for  separating  two  simple  polygons  by 
a  single  translation;  this  is  the  planar  case  of  problem  1  above.  Similarly,  Arkin, 
Connelly  and  Mitchell  [1]  address  the  planar  version  of  problem  2.  They  use  the 
concept  of  monotone  paths  among  polygonal  obstacles  to  identify  a  removable 
subassembly  of  simple  polygons  in  the  plane.  The  methods  in  [1]  do  not  extend 
directly  to  the  three-dimensional  case.  However,  Mitchell  has  independently 
shown  that  directions  for  partitioning  an  assembly  can  be  found  in  polynomial 
time  [10].  In  this  paper  we  give  an  algorithm  to  find  the  set  of  translations  sep¬ 
arating  two  polyhedra  in  0(n‘‘)  time  and  show  that  this  is  optimal.  A  method 
for  finding  complete  assembly  sequences  based  on  this  algorithm  is  analyzed 
theoretically  as  well  as  empirically  and  improvements  for  applications  are  de¬ 
scribed.  The  analytic  bounds  derived  consider  the  number  of  parts,  the  number 
of  vertices  in  the  representation  of  the  parts,  and  the  size  of  the  coordinates  in 
the  input  description. 


2  Separating  Two  Polyhedra 

In  this  section  a  method  for  finding  the  set  of  translations  separating  two  poly¬ 
hedra  will  be  derived.  Let  P  and  Q  be  closed  and  disjoint  polyhedra  in  given 
spatial  placement,  and  let  P  and  Q  be  represented  as  unions  of  at  most  n 
tetrahedra,  i.e.  P  =  lj;_(  Tj  and  Q  =  (J*  -lUi  where  r,s  <  n.  A  translation 
separating  P  from  Q  is  a  vector  d  such  that  P  can  be  translated  to  infinity  in 
direction  d  without  intersecting  Q,  i.e.  p-(-<d  is  not  in  Q  for  each  point  p  in  P 
and  each  t  in  [0,  oo). 

The  set  of  all  translations  of  P  can  be  represented  by  the  points  on  the 
unit  sphere  in  three-dimensional  space.  For  each  pair  of  tetrahedra  Ti,Uj> 
the  configuration  obstacle  C{Ti,Uj)  is  the  set  of  placements  of  Ti  such  that  Ti 
intersects  Uj  [9].  The  set  of  translations  along  which  Ti  collides  with  Uj  is  the 
projection  of  C(Ti,  Uj)  on  the  unit  sphere;  let  Rij  denote  this  region. 

The  regions  Rij  are  bounded  by  segments  of  great  circles  on  .  The  set  of 
great  circles  for  all  Rij  determines  open  regions  on  called  faces.  Each  face  is 
a  maximal  connected  component  on  the  sphere  not  intersecting  any  great  circle 
in  this  set.  The  faces  are  regular  in  the  following  sense:  the  pairs  of  tetrahedra 
Ti ,  Uj  from  P  and  Q  that  collide  in  direction  d  are  constant  for  all  translations 
d  in  a  face  /.  Let  p(/)  be  the  number  of  pairs  of  tetrahedra  that  collide  along 
translations  in  face  /.  If  the  segment  e  of  a  great  circle  lies  between  two  faces 
/  and  /',  we  have  the  following  crossing  rules  for  p(/)  and  p(/'): 

•  If  e  belongs  to  the  boundary  of  a  region  Rij  and  /'  is  on  the  interior  side 
of  e,  then  p(f')  =  p(/)  -1-  1. 
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•  If  /'  is  on  the  exterior  side  of  region  boundary  e,  then  p(/')  =  p(/)  —  1. 

•  If  e  is  not  on  the  boundary  of  any  region  then  p(/')  =  p(/). 

In  some  cases,  several  edges  may  coincide.  If  two  faces  /  and  /'  are  separated 
by  an  edge  c  bounding  several  regions,  then  p(/')  =  p(/)  —  g  +  h,  where  g  is 
the  number  of  regions  on  the  same  side  of  e  as  /  and  h  is  the  number  of  regions 
on  the  same  side  of  e  as  /'.  The  set  of  translations  separating  P  from  Q  is  the 
union  of  all  faces  /  for  which  p(f)  =  0. 

In  the  algorithm  below,  translations  are  represented  as  points  on  two  parallel 
planes  instead  of  points  on  the  unit  sphere.  Configuration  obstacles  C(Ti,Uj) 
are  projected  to  the  planes  z  =  I  and  z  =  —I  using  a  central  projection  from  the 
origin.  The  regions  Rij  arc  planar  regions  bounded  by  line  segments  and  rays. 
The  supporting  lines  of  these  segments  and  rays  define  an  arrangement  in  each 
plane,  represented  by  a  graph.  The  nodes  in  the  graph  represent  faces,  edges, 
and  vertices  of  the  arrangement,  and  links  connect  adjacent  elements.  Edges  on 
region  boundaries  are  oriented  with  respect  to  the  interior  of  the  corresponding 
region,  while  edges  obtained  by  extended  supporting  lines  are  marked  as  such. 

The  algorithm  to  find  all  faces  representing  valid  translations  proceeds  as 
follows.  For  each  plane  2=1  and  2  =  —  1, 

1.  For  each  pair  of  tetrahedra  7<,  t/j,  compute  the  projection  Rij  o{ C(Ti,Uj) 
on  the  plane. 

2.  Calculate  the  arrangement  of  lines  determined  by  the  boundaries  of  the 
regions  Rij ,  orienting  the  edges  as  described  above. 

3.  For  an  arbitrarily  selected  face  /o,  compute  the  number  p(/o)  of  regions 
Rij  containing  /q. 

4.  Perform  a  depth-first  traversal  of  all  the  faces  in  the  arrangement  by  step¬ 
ping  from  /o  to  neighboring  faces.  To  step  from  a  face  /  to  a  neighboring 
face  /',  calculate  p(/')  from  p(/)  using  the  crossing  rules  above.  After 
visiting  a  face,  it  is  marked  and  not  visited  again.  For  each  face  /  where 
p(/)  =  0,  output  the  face  /  and  continue. 

Since  P  and  Q  consist  of  at  most  n  tetrahedra  each,  there  are  at  most  n~ 
regions  Rij,  each  with  a  constant  number  of  edges.  Therefore  step  1  requires 
O(n^)  operations.  An  arrangement  of  m  lines  in  the  plane  can  be  computed  in 
O(m^)  time  [2,  4]  and  has  0{ni?)  cells.  Here  m  =  so  the  number  of  cells 
and  the  computing  time  for  step  2  are  O(n^).  Each  region  Rij  has  a  constant 
number  of  edges,  so  testing  the  initial  face  fo  for  inclusion  in  all  regions  requires 
O(n^)  operations.  Finally,  the  depth-first  search  steps  over  each  edge  at  most 
twice,  each  step  taking  constant  time.  Since  the  number  of  arrangement  edges 
is  O(n^),  the  computing  time  for  step  4  is  0(n‘’). 
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Figure  3:  Polygons  from  Pollack  et.  al.  [12] 


The  optimality  of  this  algorithm  directly  follows  from  an  example  given  by 
Pollack,  Sharir  and  Sifrony  [12].  The  example  in  [12]  concerns  two  polygons  P 
and  Q  with  r  and  s  edges  respectively;  the  number  of  connected  components 
in  the  complement  of  the  configuration  obstacle  corresponding  to  P  and  Q  is 
proportional  to  r^s'  (figure  3).  In  our  case  the  polygons  P  and  Q  are  regarded 
as  polyhedra  of  zero  volume,  and  r,  s  =  n;  the  following  holds  equally  if  P  and 
Q  are  polyhedra  with  sufficiently  small  thickness.  We  place  P  in  a  plane  p 
and  Q  in  a  plane  parallel  to  p,  but  distinct  from  p.  Then  the  plane  containing 
the  configuration  obstacle  of  P  with  respect  to  Q  does  not  contain  the  origin, 
so  the  projection  of  the  configuration  obstacle  of  P  with  respect  to  Q  on  the 
sphere  partitions  into  connected  components.  Therefore  the  set  of 

translations  separating  P  from  Q  consists  of  connected  components. 

Instead  of  decomposing  polyhedra  into  tetrahedra  in  the  above  method,  we 
can  decompose  the  faces  of  the  polyhedra  into  triangles,  and  find  the  config¬ 
uration  obstacles  for  pairs  of  triangles.  In  fact,  such  triangulations  are  often 
computed  in  geometric  modeling  systems.  The  faces  of  a  part  with  n  vertices 
can  be  triangulated  in  O(nlogn)  steps  into  0(n)  triangles  [5],  so  that  the  to¬ 
tal  number  of  steps  in  the  above  algorithm  remains  O(n^).  We  now  have  the 
following  lemma: 

Lemma  1  The  set  of  all  iranslaltons  separaiing  P  from  Q,  where  P  and  Q  are 
polyhedra  each  with  n  vertices  can  be  found  in  O(n^)  steps  and  this  is  optimal. 

To  find  separating  translations  in  which  parts  touch  each  other,  open  poly¬ 
hedra  can  be  considered  in  the  above  algorithm.  In  this  case  the  edge  and 
vertex  cells  of  an  arrangement  may  correspond  to  valid  directions.  A  count  p(c) 


6 


of  colliding  tetrahedra  is  associated  with  every  face,  edge,  or  vertex  cell  c  in 
the  arrangement.  The  arrangement  can  be  computed  and  the  cells  traversed  in 
O(n^)  steps. 


3  Partitioning  an  Assembly 

Problem  2  concerns  an  assembly  A  of  polyhedra  Pi,. Pt.  The  method  of  the 
previous  section  can  be  extended  to  find  a  translation  d  and  a  proper  subassem¬ 
bly  5  of  such  that  d  separates  5  from  the  remaining  parts  \  5. 

As  above,  each  pair  of  tetrahedra  Ti  and  Uj  from  different  polyhedra  define 
a  region  Rij  of  the  unit  sphere  S^.  A  directed  graph  G(f)  with  weighted  arcs  is 
associated  with  each  face  /  in  the  corresponding  arrangement  on  S^.  The  nodes 
of  each  graph  represent  the  objects  Pi, . . .,  P*.  The  weight  of  an  arc  from  P,  to 
Pj  in  G{f)  is  the  number  of  pairwise  intersections  of  tetrahedra  from  Pi  and  Pj 
during  any  translation  d  in  /.  Arcs  with  weight  zero  are  removed  from  G(f). 

The  graphs  G(f)  and  G(f')  for  neighboring  faces  /  and  /'  sharing  an  edge 
c  are  related  by  the  following  crossing  rules; 

•  If  e  is  a  boundary  segment  of  a  projected  configuration  obstacle  from 
tetrahedra  in  Pj  and  Pj,  and  /'  is  in  the  interior  of  the  region,  then  the 
weight  of  the  arc  from  Pi  to  Pj  is  one  greater  in  G(f')  than  in  G(f). 

•  If  e  is  a  boundary  segment  and  /'  is  outside  of  this  region,  the  weight  of 
the  arc  from  P,  to  Pj  is  one  less  in  G{f')  than  in  G(/). 

•  If  e  is  the  extension  of  a  boundary  segment,  G(f')  =  G{f). 

Similar  to  the  case  of  two  polyhedra,  if  several  edges  coincide  then  G(/')  differs 
from  Gif)  by  the  sum  of  the  changes  for  the  coinciding  edges. 

A  proper  subset  5  of  A  can  be  removed  along  a  direction  d  in  a  face  /  if 
and  only  if  there  are  no  arcs  in  G(/)  from  nodes  in  5  to  nodes  in  A  \  5.  A 
node  Pj  is  a  successor  of  Pj  in  the  graph  G(f)  if  »  =  j  or  there  is  a  path  in 
G(f)  from  Pi  to  Pj.  The  predecessors  of  a  node  are  defined  similarly.  If  the  set 
of  successors  of  every  node  in  G(/)  is  the  entire  set  of  graph  nodes,  then  there 
is  no  subassembly  that  can  be  removed  using  directions  in  face  /.  However,  it 
suffices  to  compute  the  sets  of  successors  and  predecessors  of  a  single  arbitrary 
node  Pi: 

•  If  the  set  of  successors  and  the  set  of  predecessors  of  Pi  are  both  equal 
to  A,  then  there  is  no  proper  subassembly  of  A  that  can  be  removed  in  a 
direction  in  /.  This  follows  from  the  transitivity  of  the  successor  relation. 

•  If  the  set  of  successors  Si  of  Pi  is  a  proper  subset  of  A,  then  Si  is  a 
removable  subassembly  of  A. 
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•  If  the  set  of  predecessors  S2  of  P\  is  a  proper  subset  of  A,  then  by  definition 
no  arcs  connect  nodes  in  \  52  to  nodes  in  52-  Therefore  \  52  is  a 
removable  subassembly  of  A. 

To  find  a  removable  subassembly  of  A,  we  again  project  the  configuration 
obstacles  C(Ti,  Uj)  onto  two  planes  2=1  and  2  =  —1.  However,  if  a  translation 
d  separates  a  subassembly  5  from  /I  \  5,  then  — d  separates  the  subassembly 
A\S  from  5.  Thus  it  suffices  to  search  only  one  planar  arrangement.  This  gives 
rise  to  the  following  algorithm  for  finding  a  removable  subassembly. 

1.  Calculate  the  arrangement  of  regions  Rij  on  the  plane  2=1. 

2.  Compute  the  graph  G(fo)  for  an  initial  face  fo  of  the  arrangement. 

3.  Perform  a  depth-first  traversal  over  the  arrangement,  computing  G{f}  for 
each  new  face  /.  If  in  any  graph  G(f)  the  successors  or  predecessors  of 
Pi  are  a  proper  subset  5  of  A,  output  5  and  a  translation  d  in  /. 

The  arrangement  of  projected  regions  can  be  calculated  in  0(n‘*)  steps,  and 
the  initial  graph  G(/o)  can  be  found  in  O(n^)  steps.  Finding  the  set  of  successors 
or  predecessors  of  a  node  in  one  graph  requires  O(k^)  steps.  There  are  Oin"*) 
faces  in  the  arrangement,  so  traversing  them  all  requires  0{k^n^)  operations. 
VVe  now  have  the  following  lemma: 

Lemma  2  Let  A=  {Pi,.-.,Pifc}  be  a  set  of  k  polyhedra  with  a  total  of  n  ver¬ 
tices.  it  can  be  decided  in  0{k^n^)  steps  whether  there  is  a  proper  subassembly 
of  A  that  can  be  translated  to  infinity  without  intersecting  the  remaining  parts. 
An  appropriate  subassembly  and  direction  can  be  computed  in  the  same  number 
of  steps. 

As  an  example,  consider  the  simple  configuration  of  four  cubes  aligned  along 
the  x-axis  in  figure  2.  The  corresponding  planar  arrangement  consists  of  12 
polygons  in  each  of  the  planes  2=1  and  2  =  —  1;  several  of  these  polygons 
coincide.  Figure  4  shows  the  plane  2=1.  The  projected  configuration  obstacle 
corresponding  to  cubes  Pj  and  P4  is  the  region  /?(1,4)  and  is  bounded  by  a  line 
segment  and  two  rays. 

Figure  5a  shows  the  graph  G(/?(l,4)).  P(l,4)  is  contained  in  P(l,2)  and 
P(l,3),  so  there  are  arcs  in  the  graph  from  node  1  to  nodes  2,  3,  and  4,  each 
of  weight  1.  /2(1,4)  is  contained  in  /Z{2,4),  /?(2,3),  and  /2(3,4).  Since  node  4 
has  no  successors,  it  is  a  removable  subassembly  for  translations  in  P(l,4).  If 
cubes  P2  and  P4  represent  a  single  part  P24,  the  graph  in  figure  5b  results. 
Nodes  24  and  3  form  a  strongly  connected  component,  so  cubes  2,  3,  and  4 
must  be  removed  simultaneously  for  translations  in  P(l,4). 
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Figure  5:  Graphs  for  region  /?(1,4)  where  (a)  Fi, . . . ,  P4  can  be  moved  indepen¬ 
dently  (b)  F2  and  F4  must  be  moved  simultaneously 

4  Finding  Assembly  Sequences 

The  above  method  can  be  used  to  decide  whether  there  is  a  complete  assembly 
sequence  for  an  object  with  polyhedral  parts.  Here  each  subassembly  can  only 
be  removed  by  a  single  translation,  but  each  translation  in  the  sequence  can 
involve  one  or  more  pa'-ts. 

Lemma  3  I1  can  be  decided  in  0(k^n'*)  steps  whether  polyhedra  Pi, . . . ,  Pk  can 
be  separated  completely  using  motions  where  each  subassembly  is  removed  from 
the  remaining  objects  by  a  single  translation. 

Proof:  Assume  that  A  =  {Pi, . . . ,  Pk)  can  be  disassembled.  Applying  the 
method  of  the  previous  lemma  to  A  gives  two  subassemblies,  each  consisting  of 
one  or  more  parts.  Each  application  increases  the  number  of  subeissembiies  by 
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one,  and  the  final  number  of  subassemblies  is  k.  Therefore  the  above  method  is 
applied  k  —  1  times. □ 

Finally,  let  d  be  a  bound  for  the  number  of  binary  digits  used  to  represent 
the  coordinates  of  vertices  in  the  input  assembly.  The  size  of  all  intermediate 
values  occurring  in  the  computation  is  bounded  by  0(d).  Thus  using  standard 
algorithms  for  rational  arithmetic  we  obtain  0{k^n^d?)  as  a  time  bound  for 
finding  an  ^lssembly  sequence  with  the  above  algorithm.  Here  all  computations 
can  be  perform<'d  without  loss  of  accuracy. 


5  Experimental  Evaluation 

The  above  methods  were  iriiplemented  in  C  on  a  DECstation  5000  using  floating¬ 
point  arithmetic,  with  a  number  of  modifications  giving  practical  improvements. 
The  program  was  tested  on  randomly  generated  assemblies  to  estimate  its  char¬ 
acteristics  in  the  average  case  and  find  practical  bounds  on  its  application. 

5.1  Implementation 

A  drawback  of  the  algorithm  above  is  the  storage  requirement:  the  arrangement 
may  take  O(n^)  space  to  store,  which  is  impractical  for  complicated  assem¬ 
blies.  Furthei  more,  the  number  of  cells  is  increased  dramatically  by  computing 
the  arrangement  of  the  supporting  lines  instead  of  just  the  boundary  segments 
themselves.  The  topological  sweep-line  algorithm  in  [2,  4]  sweeps  over  an  ar¬ 
rangement  of  m  lines  in  0(m)  space  and  optimal  O(m^)  time,  but  cannot  be 
extended  directly  to  the  case  of  line  segments  instead  of  lines. 

Our  implementation  addresses  these  problems  by  performing  a  vertical  line 
sweep  [8,  13]  over  the  arrangement  of  0{n^)  line  segments.  This  algorithm  only 
stores  O(n^)  of  the  cells  of  the  arrangement  at  one  time,  and  has  running  time 
of  0((n*  +  /)logn),  where  I  =  0(n‘‘)  is  the  number  of  intersections  between 
segments. 

An  imaginary  vertical  line  passes  over  the  arrangement.  The  cells  cut  by 
the  sweep  line  in  its  current  position  are  kept  in  a  sorted  list;  the  initial  list  is 
found  by  sorting  the  lines  by  slope.  Start  points  and  end  points  of  segments  and 
intersections  between  two  segments  are  events,  kept  in  a  priority  queue  sorted 
by  sr-value.  As  the  sweep-line  moves  from  left  to  right,  events  are  processed 
and  the  list  of  cut  cells  is  changed  accordingly.  Each  event  can  be  processed 
in  O(logm)  time,  so  the  total  running  time  is  0((m  -f  /)logm),  where  /  is  the 
number  of  intersection  events.  In  our  case  m  =  Thus  the  arrangement 
calculation  requires  0((n^  -h  /)logn)  steps  where  I  =  0(n‘*). 

The  vertical  sweep-line  algorithm  maintains  the  graph  G(f)  for  each  face  cut 
by  the  vertical  line.  The  graphs  for  faces  intersecting  the  initial  sweep-line  are 
propagated  down  from  an  initial  face  at  the  top  of  the  sweep-line.  To  process  an 
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Sweep  Line 


Figure  6:  An  intersection  event  in  the  sweep-line  algorithm 


event,  the  graph  for  a  new  face  is  calculated  by  stepping  from  the  face  above  it 
in  the  vertical  line,  as  described  in  section  3.  Thus  the  graphs  for  all  faces  in  the 
arrangement  are  calculated  and  checked  without  keeping  the  whole  arrangement 
in  memory.  The  total  computing  time  for  finding  an  appropriate  subassembly 
using  the  modified  algorithm  is  0{k^n^  -1-  n^logn). 

Figure  6  illustrates  the  processing  of  an  intersection  event.  The  interior  of 
region  /?(2,3)  is  below  edge  e\,  and  edge  62  is  the  lower  boundary  of  region 
/?(1,3).  The  graphs  for  faces  /i,  /2,  and  /a  have  already  been  computed;  all  the 
graph  links  have  weight  one.  When  the  sweep  line  processes  the  intersection  of 
ei  and  62  at  point  p,  the  face  is  entered.  Edge  62  is  between  /i  and  in 
the  new  sweep  line,  so  is  computed  by  stepping  over  €2  from  G{f\).  The 

interior  of  R(l,d)  is  above  €2,  so  Gif^)  is  obtained  from  G{f\)  by  deleting  the 
link  from  node  1  to  node  3.  Nodes  1  and  2  form  astongly  connected  component 
of  G(fn),  so  the  corresponding  parts  are  a  removable  subassembly. 

The  implementation  generates  configuration  obstacles  for  isothetic  three- 
dimensional  solids  (rectangloids  or  boxes)  instead  of  tetrahedra;  however,  the 
arrangement  computation  applies  to  the  general  ceise. 

5.2  Experiments 

To  evaluate  the  practical  computing  bounds  on  the  implementation,  n  random 
di.sjoint  boxes  were  generated  and  linked  together  to  form  k  complex  objects  for 
different  values  of  n  and  k.  Removable  subassemblies  were  identified  for  these 
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Figure  7:  An  assembly  of  four  random  boxes 


n 

k 

t 

■ngmi 

^mm 

^max 

^max 

4 

0.1 

0.2 

16210 

8 

D 

0.5 

0.6 

19644 

16 

n 

5.1 

4.7 

5.4 

30908 

32 

8 

48.7 

39.4 

48.7 

44163 

64 

8 

283.6 

0.016 

281.0 

284.2 

61816 

128 

16 

0.0042 

1120.8 

1243.7 

88264 

Table  1:  Computing  times  for  partitioning  composite  objects  consisting  of  iso- 
thetic  rectangular  solids  (units;  seconds  of  CPU-time  and  1024  Bytes) 


assemblies  using  the  described  implementation  of  the  above  method.  Figure  7 
shows  a  random  configuration  of  four  boxes. 

Table  1  shows  the  computing  times  and  storage  requirements  observed.  For 
each  value  of  n  and  k,  32  samples  were  run  and  the  average,  minimum,  and 
maximum  running  times  recorded  (<,  and  <mar.  respectively),  along  with 
the  maximum  storage  needed  (smax)-  In  all  cases  the  entire  arrangement  and  all 
graphs  were  computed  instead  of  stopping  at  the  first  removable  subassembly 
found. 


6  Conclusions 

Several  extensions  of  the  described  methods  might  be  considered.  Unconnected 
subassemblies  usually  require  more  complicated  fixtures  and  more  difficult  ma- 
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nipulation  than  do  connected  subassemblies.  As  a  result,  connected  subassem¬ 
blies  are  often  preferred  in  manufacturing  planning.  By  analyzing  a  connection 
graph  of  the  assembly,  the  above  algorithm  can  be  extended  to  generate  only 
subassemblies  which  are  connected. 

In  practice,  an  arrangement  of  fewer  segments  would  result  from  projecting 
the  configuration  obstacles  of  complete  polyhedra.  The  projected  configuration 
obstacle  for  polyhedra  Pi  and  Pj  is  the  union  of  all  projected  configuration  ob¬ 
stacles  Rij  of  two  tetrahedra  Ti,Uj  from  Pi  and  P2-  The  configuration  obstacles 
could  also  be  found  using  more  direct  methods  [8]. 

Finally,  other  types  of  motions  could  be  considered.  For  instance,  a  sequence 
of  translations  might  be  allowed  to  separate  subassemblies,  or  spatial  screw 
displacements  could  be  considered  instead  of  translations. 
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